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In this paper we present a particle-number-conserving (PNC) functional formalism to describe 
the dynamics of a cold bosonic gas. Treating the total number of particles as a constraint, whereby 
the phase invariance of the theory becomes local in time, we study this U(l) gauge theory using 
DeWitt's "gauge invariant effective action" techniques. Our functional formulation and earlier PNC 
proposals are shown to yield equivalent results to next-to-leading order in an expansion in the inverse 
powers of the total number of particles. In this more general framework we also show that earlier 
PNC proposals can be seen as different gauge (and gauge fixing condition) choices within the same 
physical theory. 

PACS numbers: 



I. INTRODUCTION 



Experimental successes in the making of Bose-Einstein condensates (BEC) ushered in a new epoch in 
atomic and optical physics where precision control of many-atom systems becomes possible. From this 
we see the rapid development of new tools and techniques where one can construct such systems designed 
for probing deeper physical issues as well as for more practical purposes (such as quantum information 
processing QIP). With these advances more sophisticated theoretical methods need be introduced to describe 
the quantum field and nonequilibrium properties of such systems, addressing the issues related to quantum 
coherence, correlations, fluctuations and noise. The demand here is to treat the dynamics of strongly 
correlated quantum systems often in the non-Markovian regimes while respecting the full quantum coherence 
of the system, meeting the constraints from the set-up and incorporating the effect of noise and fluctuations. 

Technically speaking, because any consistent description of quantum fluctuations around a BEC [l| is 
gapless [2 Q, individual Feynman graphs contributing to quantities of interest may be large. An accu- 
rate description of the full quantum dynamics requires resummation of whole classes of graphs. A familiar 
example is the hard thermal loop techniques developed to describe hot quark-gluon plasmas 3 . These non- 
perturbative approaches are most easily implemented at the level of the so-called closed time-path effective 
action (CTPEA) 0, 0, Q, closely related to the Feynman - Vernon influence functional for quantum open 
systems @- 

The CTPEA allows for two complementary lines of development. If the focus of interest is the evolution of 
quantum fluctuations and their correlations, then the CTPEA is best combined with the so-called 2 particle 
irreducible (or Cornwall - Jackiw -Tomboulis) [t| effective action to yield dynamical equations of motion 
for the correlations themselves |l0|. This is in fact the method of choice for dealing with quantum fields 
in far from equilibrium situations [11). We should stress that the formulation of the BEC dynamics as a 
hierarchy of correlations is well known |12| ; what sets the EA approach apart is its efficiency with respect 
to the implementation of nonperturbative resummation of diagrams |l3L Il4l 1 1 5| . 

If the focus of interest is the condensate itself, then the CTPEA may be used to derive a stochastic Gross 
- Pitaievsky (GP) equation, whereby quantum expectation values of the full theory may be retrieved as 
averages over the noise. In gravitation and cosmology, a similar problem has led to the derivation of the 
so-called Einstein-Langevin equation as the centerpiece of the stochastic semiclassical gravity program ^(| . 
Different approaches have been proposed ^3 to tackle the interaction of BEC with noncondensate atoms. 
While these approaches can be shown to be equivalent in the appropriate limits, derivation from the CTPEA 
enables one to keep track of the non-Markovian nature of dissipation and (colored) noise in a self-consistent 
manner. This formulation may also be applied to the treatment of decoherence of BEC in microtraps |l8j . a 
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problem similar to environment-induced decoherence in the context of structure formation via the influence 
functional and CTP techniques [T^,l20j . 

Functional methods have been introduced for the description of BEC dynamics (see, e.g., the lectures of 
Stoof 6J). The advantages of a field theory approach to many-body systems are well-known: economy in 
the formulation, flexibility in the use of self-consistency to improve on perturbation theory and the explicit 
enforcement of conservation laws. There is a vast literature on this subject based on the so - called "symmetry 
breaking" approach: the theory has a global U(l) symmetry whose breakdown signals condensation. In this 
approach particle number is conserved but not fixed: there are quantum fluctuations in both the number 
of particles in the condensate and the total number of particles in the system, only the mean numbers are 
fixed. 

Though convenient this approach contradicts with the actual experimental situation, where the total 
number of particles is fixed or, if there is loss of atoms, this needs be accounted for in real physical terms. 
Moreover, in typical experiments the total number of particles (several thousands) is not large enough to 
make the difference negligible. With precision experimentation, especially for QIP purposes, one should seek 
a more accurate method. This is the primary motivation for the present work. 

To the best of our knowledge, the first proposal for handling a BEC with a fixed number of particles 
has been put forward by Arnowitt and Girardea u l21l . l eading to the development of the so-called particle 
number conserving (PNC) method p3 . 123I l24l I2M l26t l27| . The goal of this paper is to formulate a functional 
PNC approach, and to show the equivalence of the functional and those earlier proposals to next to leading 
order (NLO) in an expansion in inverse powers of the total number of particles. Let us first see what the 
problem entails with some technical help. 



In recent years, field theory techniques have been systematically applied to the dynamics of BECs [ij, |28| . 
The most common approach is based on second - quantization formalism for a many-atom system. The 
field operator [x, t) which removes an atom at the location x at times t obeys the canonical commutation 
relations: 



A. Field theory approach to BECs 



[* {x,t),*{y,t)]=0 



(1) 



[* (x,t) (y,t)] = S(x-y) 



(2) 



The dynamics of this field is given by the Heisenberg equations of motion 



(3) 



where H is the Hamiltonian 




(4) 




(5) 



at 



(G) 
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The total particle number operator commutes with the Hamiltonian and is therefore conserved. The 
theory is invariant under a global phase change of the field operator 



(7) 



B. Problems with the symmetry breaking approach 



In the symmetry breaking (SB) approach to BECs, condensation is signaled by the spontaneous breakdown 
of phase invariance Eq. (JJJ, whereby \& develops a nonzero expectation value $ (c- number). We can therefore 
employ a background field separation for \P 0, l29l l30| 



* = $ 



SB 



(8) 



where $ q (q- number) is the field operator corresponding to quantum fluctuations. Various approaches differ 
in how to handle the dynamics of these two constituents. 

For later reference, we point out that a priori ^5> q is not necessarily orthogonal to $sb; also §sb is not 
necessarily simply proportional to the condensate wave function as defined below. 

A common feature of these approaches is that the total particle number 



N = / d d x 



(9) 



is not fixed. For example, let us assume that the condensate is confined within a homogeneous box of volume 
V, condensation occurring in the lowest (translation invariant) mode. Let ag be the operator that destroys 

an atom in the k mode. Then we may approximate (see more careful discussion below) 



(10) 



k^O 



Even if we treat ^f q as a linear perturbation on the condensate, the Hamiltonian is not diagonal on 
the ar. To diagonalize it, we must introduce phonon destruction operators bj: and perform a Bogolubov 
transformation 



a %. = a kb z + /3fc& + t 



(11) 



At zero temperature, the state is the phonon vacuum, |0) = for all k ^ 0. We find 



where 



but 



(N) = j d d x = V [|$sb| 2 + 



V ^ \ k *>/ V 



(12) 



(13) 



(N 2 ) = V 2 



$ss| 2 ) 2 + |*ss| 2 Un+ i ) + $f B m + $ 2 SB m* 



(14) 
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where 

The Bogoliubov coefheients a,k and (3k cannot be equal, because the canonical (bosonic) commutation 
relations imply — \[3k\ = 1, and so also fh ^ h. We conclude that necessarily (N 2 ) ^ (N) in the 

symmetry breaking approach, signaling the presence of particle number fluctuations. 



C. Particle number conserving functional approach: Global versus local gauge 

The Heisenberg equation of motion Eq. is also the classical equation of motion derived from the action 

s = J d d+1 xim*^ - J dtn (i6) 

The quantum theory of the BEC may be regarded as the quantization of the nonrelativistic classical field 
theory defined by the action functional Eq. (|16|l . where the canonical variables are iff (x, t) and its conjugate 
momentum ifriff*. This theory conserves particle number Eq. 0, and we are interested in the case in which 
particle number takes on a definite value N. We may reinforce this point by adding a constraint on the 
theory. This is achieved by introducing a Lagrange multiplier [i q (t) , thereby writing the action as 



S = d d+ 



im*—y + hn q (t) 



T*T N 

V 



- dtn 



(17) 



The original action Eq. I|16l) is invariant under a global transformation Eq. (JJJ but the new action Eq. 
(I17fl is invariant under the local (in time) transformations (a familiar theory with local U(l) gauge symmetry 
is electromagnetism) 



* e «(t) $j $t _> e-wCt)^ + ™ (18) 

provided 9 vanishes both at the initial and final times (when 9 is infinitesimal, these are just canonical 
transformations generated by the constraint). Therefore it must be quantized using the methods developed 
for gauge theories, such as the Fadeev-Popov method (See Appendix B). 

In short, the action functionals Eq. I|16(l and l|17fl lead to the same classical theories, but to inequivalent 
quantization schemes. The assertion is that Eq. I|17fl gives a more complete description of the fundamental 
theory, as it shows explicitly the gauge dependence while the result of the quantum theory is gauge-invariant. 
As it has been shown by DeWitt and others, it is possible to introduce an action functional whose variation 
yields the equations of motion for the expectation values of the fields and preserves the gauge symmetry. 
The explicit construction of this so-called "gauge invariant effective action" (GIEA) [23] for the description 
of BEC dynamics is the subject of this paper. 



D. Specific goals and what are attained 

In this paper we shall construct the gauge-invariant effective action (GIEA) for the dynamics of BECs with 
a fixed total particle number. The variation of this GIEA yields the equations of motion for the background 
fields. 

Since physical description of results depends on the choice of gauge many existing approaches which work 
within a specific gauge produce results which cannot easily be compared with each other. In the GIEA 
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approach we can place them in the same framework and be able to compare the similarities and difference in 
the physics. For a range of existing theories with particle number conserving constraints pH. 1231 123. 125l l2r| 
we show that they are indeed within the same physical theory. Most importantly the gauge-invariant 
formalism provides the necessary foundation to carry on these approaches beyond the next to leading order 
approximation in a consistent way. It is in this more ambitious program that the power of the gauge invariant 
approach becomes decisive, not optional. 

Our goal in this paper is to construct the basic structure for such a program, putting premium emphasis on 
the gauge-invariance aspect but limiting our reach only to the "one-particle irreducible" (1PI) level j32, as 
opposed to a 2PI or $- derivable approach [9J- The extension to the <&- derivable formulation (or even further, 
to the nPI effective action 33]) is necessary to match the accuracy demanded by present d ay expe riments. 
A 2PI non-gauge theory description of the nonequilibrium dynamics of BECs is presented in [lj, fl4L Il5| ; the 
2PI approach in gauge theories at large is discussed in |34| . 

Another restriction is that we shall work with the IN-OUT formulation [33] , as opposed to the closed - time 
path (CTP) or Schwinger - Keldysh formulation |5[. Beyond the next-to- leading order a truly dynamical 
theory requires a CTP formulation which yields real and causal equations of motion (3{| . Formally the CTP 
dynamical equations may be derived from those in this paper by considering the time variable to be defined 
on a closed time path ranging from the distant past to the far future and doubling back to the distant past. 
For further discussion we refer the reader to Refs. (l3Lll4| 

This paper is organized as follows. In Section II we develop a PNC approach in canonical terms, which 
will be later used as a template for the functional approach. Our presentation is close to p^j . with some 
minor technical differences (explained in Appendix A) The essence of the paper is in Section III, where we 
construct the GIEA for the description of BECs with a fixed number of particles. We derive the equations 
of motion in a large N expansion, where N is the number of particles, and show their equivalence to those 
derived by canonical methods in Section II, up to NLO. (Observe that in the literature there is also a "large 
A/"" expansion, where Af is the number of fields (lHI^ . In our case, J\f = 2 from * and These expansions 
should not be confused with each other.) 

In Section IV we conclude with some brief remarks. App endix A compares our presentation of the PNC 
formalism with those of Girardeau and Arnowitt pll l22j |. Gardiner [23] and Castin and Dum |2^[. and 
Appendix B gives a review of gauge field theory quantization and the theory of the GIEA . None of the 
material in both Appendices is new, they are added for completeness and easy reference. 

II. THE PARTICLE NUMBER CONSERVING FORMALISM 

The symmetry-breaking approach described above has the disturbing feature that, strictly speaking, sym- 
metry breaking only occurs in the thermodynamic limit. We have therefore a formalism that assumes the 
number of particles is essentially infinite. Most actual experiments deal with situations where particle number 
is bounded and actually not so large (from a few hundred to a few thousand atoms). Under this circumstance 
a condensate as described above simply cannot happen. 

In this Section we shall describe an alternative formulation which is designed to deal with gases at a fixed 
particle number. We shall call this formulation the particle number conserving formalism, PNC for short. It 
was first propossed by Arnowitt and Girardeau |2l| , and it has been extended and improved by m any other 
people over the last forty years p2L |23| . Our own presentation follows Castin and Dum pi. l26l [27| , with 
some minor technical differences which are discussed in Appendix A. 

Let us begin by discussing how is it possible to speak of a BEC in a situation where there is no symmetry 
breaking. 

A. The one-body density matrix and long range coherence 

We consider as above a second - quantized Bose field >]/. The state of the many-body system is an eigenstate 
of total particle number operator Eq. @- There is no particle exchange with the environment. 
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In this case of a finite system, there is no symmetry breaking. The symmetry broken state is essentially 
a coherent state and thus a coherent superposition of states with arbitrarily large total particle number. 
Nevertheless, there are situations where there is long range coherence across the system, thus capturing the 
essential feature of the condensed states. Sometimes these situations are referred to as quasi-condensates, 
but we shall not make this distinction, and call them BECs as their symmetry - broken siblings. 

To characterize the BEC state, let us introduce the one-body density matrix [3(| 

a(x,y,t) = (&(x,t)*(y,t)) (19) 

Long range coherence appears when a fails to decay as x and y are taken apart. 
Observe that a is Hcrmitian and nonnegative, in the sense that for any function / 

d d xd d yf* (x)a(x,y,t)f{y)hO (20) 
Therefore it admits a basis of eigenfunctions 

d d x a (x, y, t) <\> a (y, t) = n a <p a (x, t) (21) 
where the eigenvalues n a are real and nonnegative. We assume the 4> a are normalized 

(4>a, 4>p) = Sap (22) 

(f,9) = J d d xf*g (23) 

and complete 

t)<f> a (y, t)=6(x-y) (24) 

a 

The field operator may be expanded in this basis 

*(a;,t) =J2 a « (*> *) ( 25 ) 

ot 

The Bose commutation relations imply 

a a (t) , afp (t) = S al 3 (26) 

The a a (t) are operators which, at time t, destroy a particle in the one-particle state a whose wavefunction 
is 4> a (x, t) . From the definition of a we find 

(al(t)ap(t)) =n a (t)S a p (27) 

Therefore the eigenvalues n a (t) are the mean number of particles in the one-body state a at time t. We also 
have the strong identity 
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^ = 5>« W a « (*) 



(28) 



Condensation occurs when one of the n a , say a = 0, becomes comparable with N itself. Then we have, 
for large separations 



<* {x, y, t) ~ n (f>o (x, t) (j>o (y, t) 



(29) 



which displays long range coherence, as expected. Here <j>o (x, t) is the condensate wave function. We must 
stress that this is the fundamental definition; cpo(x,t) is not necessarily identical to the mean field $>sb 
introduced in the symmetry - breaking approach (for one, <&sb is not normalized). 

B. The dynamics of the condensate wave function 

We shall now discuss the dynamics of the condensate wave function (j) (x, t) and the condensate occupation 
number Nq (we switch to a capital N to emphasize its macroscopic character). We envisage a situation in 
which N is finite but large, and will seek equations of motion as an expansion in inverse powers of N. In 
preparation for this, it is convenient to scale the interaction term, writing U = u/N. 

The idea is to exploit the fact that ao is "large" and the other destruction operators are not. Let us define 
new operators A Q , a^O, through 



The X a do not commute with a. 



Vn' 



(30) 



« 



+ x a = o 



but they commute with N. Let us define the non-condensate held 



13^0 



so the Heisenberg held operator becomes 



* (x, t) = a (t) 



we have the identities 



(31) 



(32) 



(33) 



(0o, A) =0 



(34) 



aJaoA (x, t))=0 



(35) 



The condensate occupation number is 
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N = N 



(36) 



We seek a solution of the Heisenberg equations of motion Eq. © of the form 



(37) 



a = ^Ne-^ 



7i A7 
7 JVV2 JV 



(38) 



(39) 



A = A 



AA 

NV2 



(40) 



where 7 is a constant q- number and 7^7 = 77+ = 1. Expanding the canonical commutation relations we get 



7,7i 



[71,7+] = [A ,7 + ] =0 



(41) 



71,71 



[7,A 7 t] + [A 7)7 t] = 1 



(42) 



The expansion of the Heisenberg equations in inverse powers of N yields at leading order the Gross 
Pitaievsky equation (GPE) for <I> 



= i^- - H$ - u$*$ 2 + u$ 
at 



(43) 



with a possibly time-dependent chemical potential 



do 

dt 



(44) 



The actual value of /j, is derived by consistency with the normalization of $. The next-to-leading order 
(NLO) terms read 



= i7 t^$ + i^l _ HA Q u \( 2<F$ - d) A + $2 A t 

;)( L V U/ 



dt dt 



(45) 



Projecting over $ we get 



= i^^p- - ($,<, 
dt V ' 



$*$A + $ 2 A?, 



(46) 
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and finally 



dA 

= i—^ — HAq — u ($*$) A + fiA - uQ 
at 



(47) 



where 



Q[f] = /-*(*,/) (48) 
Observe that the equation for 71 is consistent with the requirement that Nq/N = 1 — (N^ 1 ) , since 



d 
df 



7 t 7i + 717 



= 



(49) 



so we may impose the condition 7^71 + 7^7 = 0. 

We see that these equations do not determine 7. Since the one body density matrix is independent of it, 
7 may be chosen freely. The simplest choice is 7 = 1, in which case Eq. 1421) reduces to 



7i,7i 



(7 = 1) 



(50) 



With this choice, the requirement that 71 + 7} = must be understood as a restriction on the allowed 
states, rather than as a strong identity, since it would conflict with Eq. JSUJl. We shall discuss this issue in 
more detail in the next section. 



III. GAUGE INVARIANT EFFECTIVE ACTION FOR BOSE- EINSTEIN CONDENSATES 



Consider the vacuum persistence amplitude for a theory with classical action S as in Eq. l|16fl . expressed 
as a path integral over paths with total particle number N 



Z Q = J e lS ' h IT ^ (/ d<ix 
Exponentiate the delta functions 



N 



(51) 



Z = I D^Dn q e lS » /h 



(52) 



s„ = s + n J d d+1 x n q (t) 



N 
V 



(53) 



Observe that now the path integral is redundant, since we may transform the fields as in Eq. I|18[l . We may 



fix the redundancy by factoring out the gauge group. Choose some function fe = f 
dfs/dO 0. Then 



such that 



1= / ^d8 6(f e -c) 



(54) 
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Inserting this into the vacuum persistence amplitude and average over c with a weight e lc / 2<T we get 

'Sfo 



Z Q = 6 / DVD[j, q e iS ^-' /h Det 



9=0 



where 

is the volume of the gauge group we wish to factor out. 



= s + h J dd+lx N (*) 



r 



2^ 



The determinant is expressed as a path integral over Grassmann fields £, ?7 



Det 



\Sfe~ 




56 _ 


-j 

9=0 J 



D(Dj] e 



B 4S-n 



(55) 



(56) 



(57) 



(58) 



Thus, we have a theory of fields 'J, ^\ fi q , (, and r\. Suppose we wish now to compute the one loop 
effective action for this theory. We redefine \& = $ + ^ q as in Eq. (JSJ). Similarly redefine ^ = $* + ^J, and 
/i g = /ig + A. Doing so introduces an undesirable cross term 



h J d d+1 x\{t) + 



To eliminate this term, we choose 



(59) 



/ = -a Ql q - H q ) + I (tx 

The ghost action now becomes 

d 



Sghost — i j dt C, 
and it decouples at one-loop order. 



+ - 2 |$| 



o-+ j d d x 



(60) 



(61) 



A. Inverse particle number as a small parameter 

The loop expansion is a formal development in which Feynman graphs are classified according to their 
topology. This has strong heuristic value but is hard to estimate a priori which order in perturbation theory 
must be reached to guarantee any prescribed accuracy. The inverse of the total particle number N provides a 
small parameter with a clear physical meaning, and so it is convenient to organize the perturbative expansion 
as a development in inverse powers of N. As we shall see, this may be achieved through a Stratonovich 
transformation of the interaction term. 

Let us begin by rewriting the classical action as (from now on, we assume h = 1) 
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S N = [ dtd d x (<*tu t $ _ - — * t2 * 2 - N^X 

J X 2N V J 



(62) 



with /i g as connection in the covariant derivative 



(63) 



We have assumed the nonlinear term scales with TV in a particular way (U —> u/N) to make sure the large 
N limit is well defined. 
We add a term 



(64) 



to the action, whose only effect is to multiply the generating functional by a constant factor, and integrate 
over the new auxiliary fields f2 and f^. The action now reads 



S N = J dtd d x |i*t£) t $ _ q,i H q, _ i (fit ^2 + + ^-0*0 - ATy- j 



(65) 



We add to this action the gauge fixing term as in Eq. I|57|) . The gauge fixing condition is given by Eq. i|60[l , 
where the gauge parameter is chosen to be a — Ns, and the ghost term Eq. I|61l) . 

Finally we rescale $ = \fN(j>, = \fN(j>* and H> q — VNip, ^ q = y/ffffi, whereby we get 



(66) 



s N [*, o, c, $] = NSi 6 + ^, n, ng, C, v, <f> 



(67) 



Henceforth, we shall write Si = S. 

In the rescaled theory there are only cubic interactions; all vertices scale like N, and all propagators scale 
as iV -1 . It follows that the power of N for any given graph is simply 1 — £, where I is the number of loops. 

B. Leading order equations of motion 

Introduce background fields Cl and \x q for the fields £7 = Q + uj, \i q = p, q + A. To leading order in N, the 
theory reduces to the classical one, with equations of motion 



iD t (j>- H4>-n<t>* =0 



(68) 



/ 



d d x 6*6 = 1 



(69) 



(70) 
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where D = Mr — i/i q is the covariant derivative evaluated with fi q as connection. To make contact with the 
PNC formalism above, we assume (cfr. Eqs. (|38|l and (|39(l ) 



-i@(t) 



A0 
~N~ 



(71) 



and ^„ = 0. We therefore get the GPE Eq. (g3J) for $, with /i = <i6/<it as in Eq. O above. 

C. The fluctuation fields at leading order 

To compute the first order corrections, namely, terms of order N° in the effective action, we must keep 
quadratic terms in the fluctuations, namely 



s N [*, n, n v c, r), *] = y>, /i„ c, ??, <t>] + NS quad 



(72) 



Sow = dtd x <i 



(jxjj'tp + (/)* ip' Ul 



-— I ,n { I ,r.,- 



i 



H LJ'U) 

2u 



<j)*lp + ■tp'cj) 



(73) 



The ghost and gauge action decouple to lowest order, but they must be included at higher orders. To 
eliminate the toip cross terms, we redefine 



U) = w + 2u<j)i}} 



(74) 



to get 



ld = J dtd d x LftDtij) - $Hi> - i (n^ 2 + ft 2 n") - 2u(j)*(j)^ij + y wXi 



+ j- I ,11 { I d d x 



(75) 



The w, A and the ghost fields are decoupled and play no further role at NLO. Making one last scaling and 
a phase shift 

= 7V- 1 / 2 e - ie V, (76) 
we find that to leading order the field ip is a quantized linear field obeying the Heisenberg equation of motion 



iD t ip - Hip - u® 2 ^ - 2u<&*<&ip 



d d y + (y,t) 



(77) 



and commutation relations 



[i/)(x,t),il>1 (y,t)] =5{x-y) 



(78) 
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We have assumed p> q = 0; D t = ^ — ifi is the covariant derivative with p, as connection. 

Let us split the fluctuation field into its components along the mean field (longitudinal) and orthogonal 
to it (transverse) 



ip (x, t) = iPo (t) $ (x, t) + f/> ± (x, t) 
where ip = ($,ip) and tp± = Qtp (cfr. Eq. (gSJl). The ETCCR imply 



(79) 



(80) 



= 



(81) 



1>±{x,t),il>{(3t,tj\=Q6{x-v) (82) 
We substitute Eq. I|79|) into the wave equation for the fluctuations Eq. (|T7|> ; using also the GPE Eq. (|43|) 



i— + A 
dt 



where 



(83) 



A = 

s 



(84) 



Projecting over <I> we get the equation for the longitudinal mode 



■ <#o . „ 



1p0 + tPo = Co [ipA 



where 



(85) 



(86) 



A = ($,A) 



(87) 



C [ip±] = ($,C[Vm.]) = u(#, 
Now observe that both A and Cq [if>±] are real, so writing 



with ETCCR 



■00 = C + «7 



(89) 
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[& *l\ = \ (9°) 

we get 

^ - 2A £ - Co [V±] (92) 

We realize that Eq. (|91|l is the counterpart to Eq. Q49J1. so to match the GIEA and PNC approaches we 
set £ = 0. Comparing Eq. Ij92(l to Eq. I|46|l. we see that the longitudinal fluctuation field is simply ipQ =71, 
where the choice 7 = 1 has been made. As discussed in the last section, we do not regard £ = as a strong 
condition, but rather as a restriction on allowed physical states. We shall discuss this issue in the next 
subsection. 

The equation for the orthogonal component ip±_ is now 

QC tys] = C ty x ] - $G) [ipx] = (93) 

where C and Co are defined in Eqs. lf%5|l and JHBJl, respectively. Since this is equivalent to Eq. (fTFjl . we 
identity tjj±_ = Aq. 

D. Gauge fixing dependence of the two point functions 

To complete the evaluation of the next-to-leading order (NLO) corrections to the mean fields and the 
one-body density matrix, we need the LO two point functions of the theory, that is, the expectation values 
of the products of two fluctuation fields. In this subsection we shall discuss the important issue of whether 
these two-point functions are gauge-fixing dependent, that is, whether they depend on the parameter s. Of 
course, physical observables, the one - body density matrix among them, must be gauge fixing independent. 

Observe that s has disappeared from the Heisenberg equation for tJj± Eq. (|93|l . so the two point functions 
built from it are automatically gauge fixing independent. As for the two point functions involving the 
longitudinal modes £ and 77, we translate the requirement £ = into 

(£ 2 ) = (^±> H^ + ?/0 = (94) 

while the ETCCR means 

i 



(ft "170 = 5 (95) 
Eq. implies 

j t (r) (t) ^ (y, 0) = -Co [<V>± {,t)i>± (y, f)>] (96) 

We see that there are nontrivial correlations between the longitudinal and the transverse quantum field, but 
they are gauge fixing independent. 
Finally 



15 



j t (V (*) V (f)> = \A, (t) - C [<^l (,*),, (f)>] 



A <„ (i) „ (f )) = (f ) - Co [<„ (t) ^ ( , f )>] 



These two equations show that we may write 

(V(t)v(t')) = {v(t)v(t')) + l -j\r A (r) 



(97) 
(98) 

(99) 



where (rj (t) r\ (t')) is gauge fixing independent. In particular, the gauge fixing dependence disappears in the 
coincidence limit t' = t. 

We may now proceed to evaluate the NLO mean fields. 

E. The mean fields at next to leading order 

We may write the GIEA to one loop order as 

r GJ [<p, fi, n] = ns [0, n, p. q , o, o, <j>] + r 1 p q ,n] (ioo) 



Ti = -iln 
To NLO we write the mean fields as 



DipDip* e^W^'^] 



= e -»W 










^ ^ +°( N 










Mg = 


N 



(101) 



(102) 



(103) 



(104) 



We also split <j)^ into components along the LO mean field (longitudinal) and perpendicular to it (trans- 
verse) 



The equations of motion for the NLO terms become 



*T V1 



(1) 



(105) 



+ - ({^o (t) $ (*, t) + i/>± (x, t) , £}} + 



(106) 
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v = 



(107) 



(2^ + <{^J-}>) + <^>* 2 



<^x) - — = 



(108) 



where we identified ip± — Aq and introduced v from Eq. I|37l) . 

In view of the analysis in the previous subsection, the explicitly s-dependent term in Eq. I) 106(1 vanishes. 
All other two-point functions are in the coincidence limit, and so they are gauge fixing independent. Wc 
may conclude that the NLO correction to the mean fields themselves is gauge fixing independent. This is 
the only property of these corrections we shall need in the remainder. 

F. The one-body density matrix at next to leading order 

From the analysis above, we see that in the GIEA approach we obtain the Heisenberg field operator for 
the bosonic field as 



V'o 



N 




(109) 



Observe that this is also the result of the canonical PNC approach Eq. I(33|l under the choice 7=1. Since 
it is observable, a m must be explicitly gauge independent. 
We get, using the NLO equations of motion Eq. (|107fl . 



a m (x, y) = (T {%, y) + 5a m (x, y) 



(110) 



a (x,y) = A$* (x)$(y) 



1 - 



N 



(111) 



Sa m (x,y) = $* (x) (y) + (vJV'i (tf) 

+ (x) + (x) V^o)) * (y) + (4 (a) <l>± (y) 



(112) 



To diagonalize this, we simply think of Sa m as a perturbation of a . er clearly admits <& as an eigenvector 
with eigenvalue N — v. Since ($, 5a m $) = 0, there is no correction to this eigenvalue at first order. This 
means that the condensate occupation number is Nq = N [l — v/N + o (-/V" 1 )] , as in the canonical approach, 
Eq. (|36() . The actual condensate wave function is 



= $ + ^ (y) + (4^ (y))) 



(113) 



which is not the same as the NLO mean field. Both the condensate occupation number and wave function 
are explicitly gauge independent, as expected. 

Finally we mention that, in spite of there being nontrivial mean fields, the expectation value of the bosonic 
field operators is zero, as it must be in a finite system. This comes about because the field operator transforms 
a physical state into one that does not satisfy the particle number constraint; for further discussion, see |39j |. 
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IV. DISCUSSIONS 

We have presented a functional formulation for cold atomic gases. This condition enters as a constraint in 
the path integral, thus changing the global U (1) symmetry of the model to a local (in time) one. Therefore 
the theory must be quantized by the Fadeev - Popov method. We derived a gauge-invariant effective action 
for the description of its dynamics. 

To make contact with the more familiar approaches, we draw attention to the fact that the mean field 
is not necessarily identical with the condensate wave function. However, once the mean field is found, it is 
easy to obtain the condensate wave function as an expansion in inverse powers of particle number. 

The approach presented in this paper unveils the physical equivalence of several proposals in the literature, 
which now can be seen as different gauge (and gauge fixing condition) choices within the same theory. Thus, 
if one takes the s — ► oo limit, one eliminates the fluctuations in the self-consistent, time-dependent chemical 
potential /i and the ghost field, but places no restriction on quantum fluctuations of the atomic field in the 
condensate mode. This is the road taken in [37j and |13|. In the opposite s — ► limit one forces quantum 
fluctuations in the condensate mode to vanish, but must include chemical potential and ghost fluctuations in 
higher order perturbation theory. As we can see, neither of these alternatives is the best in every conceivable 
situation, which testifies to the value of a general formalism as presented here. The other gain in this 
approach is of course to make contact with the extensive body of work on gauge theories. 

The more general nature of this new (gauge invariant) approach is the appearance of a new parameter 
s. While explicitly present in the intermediate steps it does not appear in the final result, as physical 
predictions should not depend on it (gauge-fixing independence). As there is considerable expertise in 
identifying the physically relevant, s - independent predictions of the theory, this is a lesser evil than its 
appearance suggests |38|. After all this is as old a topic as electromagnetism, told in the modern gauge 
theory language. 

As a first concrete application, in (3^ the functional PNC approach is used to study the one-body 
density matrix for a cold bosonic gas in an optical lattice, near the superfluid - insulator transition. In this 
application yet another gauge choice is made (in gauge theory parliance, a simpler "covariant" gauge as 
oppossed to the "background field" gauges used here). This reflects on the flexibility of the method to adapt 
to the challenges of a concrete environment, knowing that physical equivalence is guaranteed throughout by 
general theorems. 
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APPENDIX A: BRIEF SURVEY OF PNC PROPOSALS 

In this Appendix, we shall compare in some detail the PNC formalism presented in Section II with those 
of UM- We proceed in chronological order of publication. 

1. The Girardeau - Arnowitt theory 

The Girardeau - Arnowitt (GA) theory [2l| stands apart from the other two proposals to be discussed 
because its goal is to characterize the BEC ground state within a variational approximation, rather than a 
dynamical theory. Only the homogeneous case is discussed, so linear momentum is a good quantum number 
for one-particle states. Finally, GA do not perform a large N expansion; they aim to describe the BEC 
regime even when depletion is large. 

The basic insight of the GA theory is that in the Bogoliubov theory the ground state is seen as a coherent 
superposition of states with different numbers of particle pairs (a particle of momentum p correlated to one 
of momentum —p) and thus different total particle numbers (see any textbook presentation, e. g. 40] ). GA 
describe the ground state as a coherent superposition of particle - hole pairs (a particle of momentum p, 
another one of momentum —p, and two particles less in the condensate), all of them having the same total 
particle number. A particle - hole is destroyed by the operator c p — a p /3g (2^, where a p destroys a particle 
of momentum p and 



00 = (N + ir 1/2 



(Al) 



where we follow the updated formulation of j^. In terms of the notation of Section II, Eqs. I|38|l. I|39[l and 
|@5| we have (we assume 7 = 1) 



Bo 



N 



No + l 



-iB(t) 



1 



-i&(t) 



7i 



1 



7i 



A7 

A 7 r 



ATi/2 N 

From Eq. (|3U|> . the particle - hole destruction operator is 



(A2) 



c P = -^=\ p No{No + l) 1/2 
v N 



No 



^N (N + 1) 



X p + (N- 1 ) 



(A3) 



In the homogeneous case $ = <po — 1/s/V. Therefore Eq. (j46(l gives 71 = constant, the Gross - Pitaievsky 
equation l|4^|) yields [l = u/V and Eq. (P7Jl becomes 







.dX p 



2m 



-A, 



V 



[Ap + At] 



(A4) 



which also follows from the GA Hamiltonian at leading order (see Eq. (14) in |22|) 
These formulae provide the translation between the GA theory and our formalism. 
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2. The Gardiner theory 

Gardiner's theory [2j| is also a systematic large N expansion of the Heisenberg equations of motion, after 
introducing an ansatz equivalent to Eq. (|33f) for the field operator. The condensate wave function is not 
introduced as an eigenfunction of the one-body density matrix, but rather kept as an unknown. The Gross 

- Pitaievsky equation Eq. (fO)) is derived as a consistency condition, because if it does not hold, then the 
time derivative of the non-condensate field acquires a term of order TV 1 / 2 , thus invalidating the large N 
expansion. This procedure is sufficient for a NLO discussion as carried out in this paper, but the condensate 

- non condensate separation becomes ambiguous at higher orders, in which case the Castin and Dum proposal 
p3 | may be preferable. 



3. The Castin and Dum theory 

As we have stated in the Introduction, our presentation of the PNC method in Section II follows the Castin 
and Dum theory [24| • The only difference is that Castin and Dum define the non-condensate operator as (in 
the notation of Section II) 

Aexfot) = j^Afot) (A5) 

This makes Eq. <|35[l simpler 



(A ex (x,t))=0 (A6) 

allowing for a clearcut separation between the condensate and noncondensate contributions to the equations 
of motion, but the expression of the Heisenberg field in terms of A ex (x, t) is more involved than in terms of 
A (x, t) . Otherwise, both formulations are fully equivalent. 

For further developments within the PNC formulation, see [2^, |2]} . 



APPENDIX B: THE GAUGE INVARIANT EFFECTIVE ACTION 



In this section, we shall review the theory of the GIEA. We begin by introducing the effective action in a 
general non-gauge theory. 



1. The effective action | 4l| 

Consider a theory of fields whose dynamics is described by an action S [vt 4 ] . The evolution of the fields 
from some initial to some final times (usually minus and plus infinity, respectively) can be represented by 
the vacuum persistence amplitude defined as the amplitude for the initial vacuum state to evolve into the 
OUT vacuum state: 



Z [0] = (Oout \0in) = e tW[0] = j e lS ^ (Bl) 

where we set Ti — 1 and adopt the DeWitt convention that a single index i denotes both space-time and 
internal indices, and repeated indices denote integration over space-time and summation over internal indices. 
More general initial and final states may be considered, but this shall be enough for our purposes. 

If we allow the fields to interact with external c- number sources Ji, then the persistence amplitude becomes 
the generating functional. Its Legendre transform is the effective action, which obeys 



20 



e <r[»] = J e i(s[* 4 ]-r.,(*'-»')) (B2) 

Observe that 

e -ir[*] /" $ < e <(s[*«]-r tl (**-»«)) = $ < (B3) 

Consider fluctuations around the background fields with ^ — <£> l + f/; 4 , and expand the action for ip l 
to second order in 



S [$ + V] ~ 5 [$] + + W (B4) 



Then 



r [$] = 5 [$] + ST (B5) 

where 

e* 41 *!*] - J Dtf eKis.i^-sr,^) (B6) 
To leading order, ST^ on the right hand side may be neglected, then the one-loop effective action is given by 

ST = - In Det [S^]' 1 (B7) 



2. The effective action for gauge theories I42L |4 

Let us consider now the case where the classical action is invariant under an infinitesimal gauge transfor- 
mations parametrized by e a 

¥ -> ¥ [e] = * 4 + (B8) 
£4< J = Ql e a (B9) 
S,iQ l a = 0. (BIO) 

We assume the gauge transformations are linear 

Qa, jk = 0, (BH) 

that they form a group 
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QiM - QpM = c lpQ^ ( B12 ) 

and are volume-preserving 

Qa,i = (B13) 
The structure constants satisfy the Jacobi identity and 

C a a p = (B14) 

Our previous definition of the persistence amplitude is inadequate because gauge equivalent configurations 
are counted as distinct. To obtain a satisfactory definition we must modify the measure of integration in 
such a way that each gauge orbit is counted only once. 

Let us introduce linear functions f a ['J 2 ] which are not gauge invariant; in particular, we request 



- fiCfy (B15) 



to have an inverse. Then 



J De a Dct [Ajg [e]]] S [f a [¥ [e]] - c a ] = 1 (B16) 
(provided the argument of the delta functions vanishes somewhere) , and we may write 

Z[0] = J De a j D^ 1 Det [A£ [** [e]]] 5 [f a [e]] - c a ] e iS M (B17) 
Since both the classical action and the volume element are gauge invariant, this is the same as 

Z [0] = J De a J D^ 1 Dct [A£ [**]] 5 [f a [¥] - c a ] e <s M (B18) 
This means the group volume may be factored out, thus defining 

Z [0] = J D¥ Det [Af 4 [*']] 8 [f a [^] - c a ] e^M (B19) 

In this expression, the correct integration measure is displayed. 

It is interesting to check explicitly that this expression is independent of the gauge fixing conditions f a . 
Suppose we replace the functions f a by new gauge fixing conditions f Q = f a + Sf a . At the same time, we 
change variables in the functional integral to new fields defined by the condition that f a [*] = f a [*] . 
This change is actually a gauge transform, since we may choose 

* J = * J - Q\e a (B20) 

where s a — [A^ 1 ]^ Sf 13 . Then the argument of the delta function as well as the exponent remain unchanged. 
Since the gauge parameters may be field dependent, the volume element is not invariant 
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= DV l Det— = DV l 
Concerning the functional determinant, we have 



(B21) 



[*] 



/,?< 



J /3 



C* 



Since = 0, 



D^Det [A| [**]] = £>**Det [A^ [¥*]] 



(B22) 
(B23) 
(B24) 
(B25) 



(B26) 



This concludes the proof of gauge independence. 

Since the persistence amplitude is independent of the c a , we may further average over them with a Gaussian 
weight exp [ic a c a /2a]. Also the functional determinant may be expressed as a path integral over Grassmann 
fields 



Det [A a p ] = J DCDrf exp{-Q a ^} 



(B27) 



The final expression looks like an ordinary persistence amplitude for a theory with a modified action 
functional 



Z[0] 



D^ 1 DCDt] p e 



(B28) 



S. 



FP 



S[tf'] 



'la 



(B29) 



We may use this action to build an effective action in the usual way. If the gauge fixing functions / are 
chosen conveniently, it is possible to make the effective action gauge invariant as well, thereby becoming the 
GIEA. 



3. The gauge invariant effective action | 3l| 



We now belabor the dependence of the effective action on the gauge fixing functions f a . We have seen 
that the persistence amplitude is independent of the choice of gauge fixing conditions. Since the effective 
action reduces to the logarithm of the persistence amplitude when the equations of motion hold (for short, 
"on shell"), the variation of the effective action with respect to the f a must be related to the first derivatives 
of the effective action with respect to the background fields. 

To make this explicit, let us consider a change in the gauge fixing condition f a — ► f a +Sf a . As we have seen, 
the volume element, action, gauge fixing condition and ghost action remain invariant if we simultaneously 
make a gauge transformation f -> <f l - Q l a e a , where e a = [A" 1 ] " 5f . But the source terms are not 
invariant, whereby 
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sr = r 



(B30) 



This formula is the key to define a gauge invariant effective action. The idea is to allow the gauge fixing 
conditions f a to depend parametrically on the background fields <I> 1 , in such a way that the total variation 
of the effective action is orthogonal to gauge trajectories 



In view of the previous formula, this means 



sr sr sf 
s¥ + Jp~s¥ 



(B31) 



{<P]e- ir J D^ 1 



DCDrf 



Sf 



& L J S& 



i[s F p[*',c-,^]-r, »(**-**)] = (B32) 



Since the are linear, we may also write 



r 4 e- 4r /" D&DCDrf Q\ [*] [A" 1 ]^ [*] 



ol [s F p[**,c°,r,' 3 ]-r, fc (* fe -<i. fc 



(B33) 

The simplest way to satisfy this condition is to choose gauge fixing conditions f a invariant under a 
simultaneous transformation of and If fields, in which case the brackets vanish identically. Observe that 
since the total derivative of the GIEA is proportional to the partial derivative (that is, with the gauge fixing 
condition kept fixed), the equations of motion derived from one or the other are equivalent. 
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